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Abstract. The scalar difference equation Xn+i = fnixmXn-i, ■ ■ ■ ,Xn-k) may 
exliibit symmetries in its form that allow for reduction of order through sub- 
^vq I stitution or a change of variables. Such form symmetries can be defined gen- 

erally using the semiconjugate relation on a group which yields a reduction 
^/^ • of order through the semiconjugate factorization of the difference equation of 

Q . order k -\- 1 into equations of lesser orders. Different classes of equations are 

(-^ i considered including separable equations and homogeneous equations of de- 

j^ I gree 1. Applications include giving a complete factorization of the linear non- 

homogeneous difference equation of order A; -|- 1 into a system of k -\- 1 first 
order linear non-homogeneous equations in which the coefficients are the eigen- 
^S] ■ values of the higher order equation. Form symmetries are also used to explain 

^ ! the complicated multistable behavior of a separable, second order exponential 

[^^ I equation. 
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, J . 1 Introduction 

c^ ' Certain difference equations have symmetries in their expressions that allow a reduc- 

tion of their orders through substitutions of new variables. For instance, consider the 
second order, scalar difference equation 

where Qn is a real function for each integer n. This equation has a symmetry in its 
form that is easy to identify when ([T]) is re-written as 



Now, setting t„ = x„ — x„_i changes Eq.([2]) to the first order equation 

tn+l = gnitn)- (3) 

The expression x„ — x„_i is an example of what we may call a form symmetry. 
Substituting a new variable t„ for this form symmetry in ([2]) gave the lower order 
equation ([3]). The form symmetry also establishes a link between the second order 
equation and the first order one, in the sense that information about each solution 
{tn} of ([3]) can then be translated into information about the corresponding solution 
of (fl]) using the equation 

n 
Xn = Xn-\ + tn = Xo + ^ tfc (4) 

k=\ 

where xq is an initial value for ([T]). Along similar lines, the non-homogeneous linear 
difference equation 

x„+i + px„ + qXn-\ = an with 1+p + q = (5) 

has at least two form symmetries. First, setting p = —1 — q and rearranging terms 
in ([5]) reveals the form symmetry x„ — Xn^i and the corresponding order-reducing 
substitution: 

Xn+l -Xn = an + q{Xn - a;„_i) => tn+l = "n + qtn- (6) 

Further, 

X„+i - qXn = an- {p + q)Xn - qXn-1 = an + Xn- qXn-1 => tn+l = "n + ^n- (7) 

Thus Xn — qXn-i is also a form symmetry of ([5]). Note that the coefficients q and 
1 of tn in (ini) and (171), respectively, are both eigenvalues of the homogeneous part of 
([5]), i.e., roots of the characteristic polynomial z^ + pz + q when p + q = —1. Later in 
this article we show that this relationship holds for all linear difference equations. 

These and many other types of known form symmetries of difference equations of 
type 

can be defined in terms of semiconjugacy; in p2] there is a basic discussion of this 
topic for real functions but we give a more general definition in this article. The idea 
of reducing order via semiconjugacy is basically simple; we find functions that are 



semi conjugate to the given functions /„ but which have fewer variables than k + 1, 
i.e., the order of ([8]). Then the semiconjugate relation simultaneously determines 
both a form symmetry and a factorization of Eq.(l8]). This factorization of difference 
equations, which is also a formal representation of the substitution process discussed 
above, yields a pair of lower order equations. This pair is made up of a factor equation 
such as ([3]) and an associated equation such as (jlj) that is derived from the form 
symmetry and relates the factor equation to the original one. The orders of the 
factor equation and the associated one always add up to the order of the original 
scalar equation. 

The aim of this article is to formalize, within the framework of semiconjugacy, 
the concept of form symmetry and its use in reduction of order by substitutions. In 
addition to unifying various ad hoc techniques, this approach also gives rise to new 
methods for analyzing higher order difference equations. Some of these new methods 
are described in this article along with examples and applications. 

Symmetries of a different kind have already been used to study higher order differ- 
ence equations or systems of first order ones. A well-known approach involves adapta- 
tion of the Lie symmetry concept from differential equations to difference equations; 
see, e.g., for a general discussion of the discrete case covering various topics such 
as reduction of order, integrability and finding explicit solutions. Also see [4], fT3], 
|15] for additional ideas and techniques. The main difference between the concepts of 
form symmetry and Lie symmetry may be summed up as follows: Form symmetries 
are sought in the difference equation itself whereas Lie symmetries exist in the solu- 
tions of the equation. The existence of either type of symmetry can yield valuable 
information about the dynamics and the solutions of the difference equation with a 
variety of applications such as reduction of order. 

2 Semiconjugate forms 

The material in this section substantially extends the notions in [22] to the more 
general group context. For related results and background material on difference 
equations see, e.g., [T], [3], [5], [8], [T2] . 

In this article, G denotes a non-trivial group. The group structure provides a 
suitable framework for our results. However, in most applications G turns out to be 
a substructure of a more complex object such as a vector space, a ring or an alge- 
bra possessing a compatible or natural metric topology. In typical studies involving 
difference equations and discrete dynamical systems, G is a group of real or complex 
numbers. We may define each mapping /„ on the ambient sturcture as long as the 



following invariance condition holds 

fn{G''^^) C G for all n. (9) 

If ([9]) holds then for each set of initial values xo,x_i, . . . ,X-k in G Eq.(IH]) recur- 
sively generates a solution or orbit {xn}'^=-i in G. Before defing form symmetries for 
([8]), it is necessary to discuss some general defintions involving systems. 

Let 1 < m < k. We say that a self map F„ = [fi^n, ■ ■ ■ , fk+i,n] of G^~^^ is semicon- 
jugate to a self map $„ of G™ if there is a function if : G'''^^ -^ G"" such that for 
every n, 

HoFn = ^nOH. (10) 

Each mapping $„ is called a semi conjugate factor of the corresponding F„. Sup- 
pose that 

H{uq, ...,Uk) = [hi{Uo, . . . , Mfc), . . . , hm{uo, ■■■, Uk)] 
$„(tl, . . . , tm) = [4'l,n{tl, • • • , tm), • • • , 4>m,n{tl, • • • , ^m)] 

where hj : G''=+^ ^ G and (f)j,n '■ G"^ -^ G are the corresponding component functions. 
Then identity (ITOl) is equivalent to the system 

hj{fl,n{Uo,- ■ -jUk),. . . , /fc+l,„(Mo, • • • ,Uk)) = (f)j,n{hl{Uo, . ■ .,Uk),. ■ ■ ,hm{Uo,. ■ ■ ,Uk)) 

j = l,2,...,m. (11) 

If the functions fj^n are given then (ITTl) is a system of functional equations whose 
solutions hj,(l)j^n give the maps H and $„. The functions $„ on G™ define a system 
with lower dimension than that defined by the functions i^„ on G'^^^. For a given 
solution {Xn} of the equation 

Xn+l = Fn{Xn), Xq E G (12) 

let y; = if (X„) for n = 0, 1, 2, . . . Then 

F„+l = ii(X„+i) = H{Fn{Xn)) = ^n{H{Xn)) = ^n{Yn) 

so that {Yn} satisfies the lower order equation 

F„+i = <l>n{Yn), Yo = H{Xo) G G"\ (13) 

The relationship between the solutions of (TT2l) and those of (TT3|) is not generally 
straightforward; however, information about the solutions of (fT3|) can shed light on the 
dynamics of (iT2l) . The case where G = R, m = 1 and i^„ = i^ is time-independent (or 
autonomous) is of some interest because in this case Eq. ( iT3l) is a first order difference 
equation on the real numbers and as such its dynamics are much better understood 
than that of the higher dimensional Eq. (ll2l) . This case is discussed in detail in [22] . 



3 Order-reducing form symmetries 

For the scalar difference equation ([8]) that is of interest in this article, each Fn is the 
associated vector map (or unfolding) of the function /„ in Eq.(IHl), i.e., 

Fn{uo,. . .,Uk) = [fn{Uo,- ■ ■ ,Mfc),Mo,. . -^Uk-l]. 

Even if each such F„ is semiconjugate to an m-dimensional map $„ as in flTO|) . 
the preceding discussion only gives the system ( fT3l) in which the maps $„ are not 
necessarily of scalar type similar to Fn- While this may be unavoidable in some cases, 
adding a few reasonable restrictions can ensure that each $„ is also of scalar type. 
To this end, define 

/ii(mo, . . . , Mfc) = Mo * h{ui, . . . , Mfc) (14) 

where h : G'' -^ G is a function to be determined and * denotes the group operation. 
This restriction on H makes sense for Eq.(IH]), which is of recursive type; i.e., x„+i 
given explicitly by functions /„. With these restrictions on H and F„ the first equation 
in (TTTj) is given by 

/„(mO, . . . , Ufc) * h{uo, . . . , Uk-l) = gn{uo * h{ui, . . . , Mfc), . . . , hm{uo, . . . , Ufc)) (15) 

where for notational convenience we have set 

9n = 4>l,n '■ G -^ G. 



Eq. fjTSj) is a functional equation in which the functions /i, hj,gn may be determined 
in terms of the given functions /„. Our aim is ultimately to extract a scalar equation 
of order m such as 

tn+l = QnV'ni ■ ■ ■ i ^n-m+l) (lo) 

from (1151) in such a way that the maps $„, will be of scalar type. The basic framework 
is already in place; let {x„} be a solution of Eq.([8]) and define 

Then the left hand side of (HM is 

which gives the initial part of the difference equation (TT6|) . In order that the right 
hand side of (TT5l) coincide with that in (fT6|) it is necessary to define 

i^jKp^ri) • • • T-^n—k) ^"n—j+1 ■^n— j'+l * '^\-^n—jy ■ ■ ■ ) •^n—k—j+l) i J Z, . . . ,771. \^ ' ) 



Since the left hand side of (TT7|) does not depend on terms Xn~k-i, ■ ■ ■ , Xn~k-j+i it 
follows that the function hj must be constant in its last few coordinates. Since h deos 
not depend on j the number of constant coordinates is found from the last function 
hm- Specifically, we have 

iT"m\^n) • • • ) ■^n—k) •^n—m+1 * i^\-^n—m) • • • ; •^n—ky •^n—k—ly • • • ; •^n—k—m+l) \^^) 

k—m+1 variables m— 1 terms with h constant 

The preceding condition leads to the necessary restrictions on h and every hj for a 
consistent derivation of flT6|) from flTSl) . so flTSj) is a consistency condition. Now from 
(fTTl) and ([IH]) we obtain for (mq, • • • , Wfe) e G'=+^ 

hjiuQ, . . .,Uk) = Uj-i * h{uj, Uj+i . . . , Uj+k-ni), j = 1, . . . , m. (19) 

We refer to H = [hi, ... , hm] as a /orm symmetry for Eq.(l8]) if the components /ij 
are defined by ( TT9l) . Since the range of if has a lower dimension than its domain, we 
say that H is an order-reducing form symmetry. 

Using the forms in (fT9l) for /ij in (fT5l) for every solution {x„} of Eq.Q we obtain 
the following pair of equations from (TTSi) . the first of which is just ( !T6|) : 

''^n+l = Qnitni ■ ■ ■ i ^n-m+l); (20a) 

X„+i = t„+l * h[Xn, . . . , Xn-k+m) ■ (20b) 

The power —1 represents group inversion in G. The first equation (I20ap may be 
called a factor of Eq.(l8]) since it is distilled from the semiconjugate factor $„. The 
second equation ( ]20b[) that links the factor to the original equation may be called 
a cof actor of Eq.(IH]). We call the system of equations ( 120|) a semiconjugate (SC) 
factorization of Eq. ([8]) . 

Note that if {t„} is a given solution of (I20al) then using this sequence in fl20bp 
produces a solution {x„} of (IE]). Conversely, if {a;„} is a solution of (IE]) then the 
sequence t„ = x„ * h{xn-i, ■ ■ ■ , Xn-k+m-i) is a solution of (I20ap with initial values 

t-j = X-j * h[X-j-i, . . . , X-j-k+m-i), J = 0, . . . , m — 1. 

Since solutions of the pair of equations (I20ap and ( I20bp coincide with the solutions 



of the scalar equation ([8]), we say that the pair (I20ap and fl20bp is equivalent to ([8]). 
The following summarizes the preceding discussions. 

Theorem 1. Let k>l, l<m<k and suppose that there are functions h : 
Qk-rn+i __^ Q ^^^ g^ . Qm __^ Q ^^^^ satisfy CQuations fTS\) and fW\} . Then with the 
order-reducing form symmetry 

H{uo, . . . ,Uk) = [Uo* h{ui, . . . , Uk+l-m), • • • , Um-l * h{Um, ■ • • , Mfc)] 



Eq. ^ is equivalent to the SC factorization consisting of the pair of equations Ii20a\) 
and Ii20b\) whose orders m and k + 1 — m respectively, add up to the order of ^. 

In this setting we say that the SC factorization fl20l) gives a type-(m, k + 1 — m) 
order reduction for Eq.(l8]), or that ([H]) is a tjpe-{m, k + 1 — m) equation. A second 
order difference equation {k = 1) can have only the order-reduction type (1,1) into 
two first order equations although the factor and cofactor equations are not uniquely 
defined. In general, a higher order difference equation may have more than one SC 
factorization. A third order equation can have two order-reduction types, namely, 
(1,2) and (2,1). Of the k possible order reduction types for an equation of order k + 1 
the two extreme ones, namely, (1, k) and {k, 1) have the extra appeal of having an 
equation of order 1 as either a factor or a cofactor. In the next two sections we discuss 
classes of higher order difference equations having one of these order-reduction types. 

We note that the SC factorization of Theorem 1 does not require the determina- 
tion of (pj^n for j > 2. For completeness, we close this section by showing that each 
coordinate function (pj^n projects into coordinate j — 1 for j > 1, thus showing that $„ 
is of scalar type, i.e., it is the unfolding of Eq. flTB]) in the same sense that F„ unfolds 
([8]) . If the maps hj are given by (TT9|) then for j > 2 (TTTj) gives 

0j>(^i(MO, • • • , Mfc), • • • , hm{Uo, ..., Mfc)) = hj{fn{Uo, . . . , Mfc), Mq, . . . , Mfc-l) 

= Uj^2 * h{Uj-i, Uj . . . , Uj+k-m-l) 

= hj_i{uo,...,Uk). 
Therefore, for each n and for every (ti, . . . , t^) G H{G^~^^) we have 

$„(tl, . . . , tm) = [gnitl, • • • , tm),ti, . . . , tm-l] 

i.e., $n|ji^(G'=+i) is of scalar type. Further, if H is defined component-wise by flT9|) then 
H{G^^^) = G^; i.e., H is onto G"^ so that $„ is of scalar type. To prove the onto 
claim, we pick arbitrary [ti, . . . , tm] G G"" and set Um-i = tm * h{um, Mm+i • • • , Wfc)~^ 
where Um = Wm+i = ■ ■ -Uk = 1 (the group identity). Then 

tm = Mm-1 * h{l,l . . .,1) 

= Um-1 * h{Um, Um+1 • • • , Mfc) 

= hmiUo,. . .,Uk) 

= hm{Uo, ..., Um-2, tm * h{l, 1 . . . , 1)"\ 1 . . . , l). 

for any choice of Kq, • • • , Mm-2 £ G- Similarly, define Um-2 = tm-i*h{um-i, Um • • • , Uk-i)~^ 



so as to get 

tm~l = Um^2 * h{Um~l, Um ■ ■ ■ , Mfe-l) 
= hm-l{Uo,...,Uk) 
= hm-l{Uo, . . . , Um-3, tm-l * h{Um-l, 1 • • • , 1)" , "m-l, 1 • • • , 1) 

for any choice of Uq, • • • , Wm-s ^ G. Continuing in this way, induction leads to selection 
of Um-i, ■ ■ ■ ,uo such that 

tj = hj{uo,. . .,Mm-i,l,- ••,!), j = 1,. . .,m 

and it is proved that H is onto G"". 

4 HDl and other type-(A;, 1) factorizations 



li m = k then the function h : G ^ G in ( 1191) is of one variable and we obtain a 
type- (A;, 1) order reduction with form symmetry 

H{uo, ...,Uk) = [uo* h{ui),ui * h{u2) ■■■, Uk-i * h{uk)] (21) 

and SC factorization 

f-n+l 9n\J^nj ^n— 1; • • • ; ^n— fc+lj 
X„+i = tn+1 * hyXn) 

where the functions Qn '■ G^ ^ G are determined by the given functions /„ in ([8]) as 
in the previous section. 

The simplest example of a non-constant h in this setting is the identity function 
h{u) = u for all u & G. An example of a type-(A;, 1) difference equation having 
this type of form symmetry over (0, oo) under ordinary multiplication is the rational 
equation 

Xn+i= ^^"-^ a,6>0. (22) 

^n^n— 1 ~r 

The term x„a;„_i in the denominator suggests multiplying fl2^ by x^ on both 
sides and substituting 

^n •^Ti'^n— 1 •^n'T'\-^n—l) 



to get the SC factorization 



Wi = - — hi = 9{in), to = Xox_i (23) 



Xn+l 



tn + b 

tn+1 ^n+1 



Xji llyXfij 



Further, Eq. (l23|) can be made hnear by the change of variables s„ = l/t„. For an 
exhaustive treatment of (!22|) based on these ideas, see [16]. Another form symmetry 
of type (I2T!) that is defined on C or M is based on h{u) = cu where c is a fixed, nonzero 
complex or real number. This type of form symmetry (with real c) has been used in 
e.g., f7] and [Hj. 

In the case where h{u) = u'^ is based on group inversion, it is possible to identify 
the class of functions /„ that have the form symmetry fl2T]) . Equation ([8]) is said to 
be homogeneous of degree 1 (HDl) if for every n = 1, 2, 3, . . . the functions /„ are 
homogeneous of degree 1 relative to the group G, i.e., 

/„(% *t,...,Uk*t) = /„(%, . . . , Mfc) * t for alH, Uj G G, i = 0,...,k, n>l. 

If G is non-commutative then this definition gives a "right version" of the HDl 
property; a "left version" can be defined analogously. We note that the two equations 
([1]) and ([5]) in the Introduction are HDl relative to the additive group of real numbers. 
For comments on homogeneous functions and their abundance on groups we refer to 
[2U] ; though stated for functions of two variables, the results in [20] easily extend to 
any number of variables. The following result shows that the HDl property charac- 
terizes the inversion-based form symmetry and yields a type-(/c, 1) order-reduction in 
every case. 

Theorem 2. Eq.^ has the inversion- based form symmetry 

H{uo,...,Uk) = [uo*u^^,...,Uk-i*Uf:^], h(t)=t'^ (24) 

if and only if fn is HDl relative to G for all n. In this case, ^ has a type-(k, 1) 
order-reduction with the SC factorization 

tn+l = /n(l, t~ , (t„ * t„_i)~ , . . . , (t„ * tn-1 * ■ ■ ■ * tn-k+l)^ ) (25a) 

X„+i = tn+1 * Xn- (25b) 



l^ote that the factor difference equation / (^5a|) has order k and its cofactor l(25b\) 
is linear non-autonomous of order one in Xn- 

9 



Proof. First, assume that ([8]) has the form symmetry (!24l) that satisfies Eg. (11 51) 
for given functions Qn, i.e., 

/n(Mo, ■ . ■,Uk)*UQ^ = gniuo * U'^'^ , . . . , Ufc-1 * u^^) . (26) 

Let t G G be arbitrary. Then for all n (!26l) implies 

/„(mO * t, . . . , Mfc * t) = 5(„((mo * t) * (mi * t)~\ . . . , (ma:-1 * t) * {Uk * t^^) * (mq * t) 
= [gniUo * Mj"\ . . . , Uk-l * U^^) * Mo] * t 
= fn{Uo,...,Uk) *t. 

It follows that /„ is HDl relative to G for all n and the first part of the theorem 
is proved. The converse is proved in a straightforward fashion; see [2T] . 

Remarks. 

1. Equation (]25bl) can be solved explicitly in terms of a solution {t„} of (I25al) as 
follows: 

Xn = YY^=otn-i * xo n = 1, 2, 3, . . . (27) 

where the multiplicative notation is used for iterations of the group opreation *. In 
additive (and commutative) notation, ( !27l) takes the form 

a;™ = xo + ZliLi^i- (28) 

2. We can quickly construct Eq. (I25ap directly from ([8]) in the HDl case by making 
the substitutions 

I ^ Xn, {tntn-l---tn-i+lY^ ^ Xn-i for Z = 1, 2, . . . , fc. (29) 

Recall that 1 represents the group identity in multiplicative notation. In additive 
notation (129!) takes the form 

U *■ ^nj ~''n ~ tn-l ■ ■ ■ — tn-i-\-l ^ Xn-i lOr Z = 1, 2, . . . , K. ("jU) 

Previous studies involving HDl equations implicitly use the idea behind Theorem 
2 above to reduce second order equations to first order ones; see e.g. [B], [101, P^ - 
Examples 1-3 next illustrate Theorem 2 and some associated concepts. 

Example 1. Consider the rational delay difference equation 

Xn+l = Xni "" ""^"^^ +bn] , Xq, X_i, . . . , X__fc > (31) 

\ Xfi—k / 

10 



where {a„}, {bn} are sequences of positive real numbers. This equation is HDl relative 
to the group (0, oo) under ordinary multiplication. Thus Theorem 2 and (!29l) give 
the SC factorization of (|3T1) as 

, _ /1^ I (^nitntn-l ■ ■ ■tn-k+2) , t \ _ j. , t 

l^n+1 — Uj I . 7 77:^^ 1- 0„ I — ttntn-k+l + 0^, 

V (,f^ntn-l . . . tn-k+l) J 

In this case, the factor equation is linear non-homogeneous with a time delay of 
k — 1 and can be solved to obtain an explicit solution of (ETl) through ([17]), if desired. 
Alternatively, we can quickly derive information about the asymptotic behavior of 
fl3T]) . For instance, if 

lim a„ = a > 0, lim bn = b > 0, a + b ^ 1 

n— >oo n— >oo 

then we conclude that all positive solutions of ^^ converge to zero if a + 6 < 1 and 
to oo if a + 6 > 1. 

Example 2. This example illustrates a situation where (IHl) and fl25ap are both HDl, 
although with respect to different groups. Consider the third order equation 

Xn+i = a;„ H " _ """^ , a 7^ 0. (32) 

■^n—l •^n—2 

Relative to the additive group G = M, this equation is HDl with the evident form 
Xn — Xn-i- Specifically, (132|) has the form symmetry 

if (m, f , w) = [u — v,v — w], h(t) = — t. 

Making the substituion t„ = a;„ — x„_i, or using (J30ll we get the SC factorization 

atl 

tn+l = - — —, to = Xq — X-i, t-i = X-i — X-2 (33) 

tn-1 

X-n+l ''n+1 ~r X^- 

This is a type-(2,l) order reduction. Note that t„ 7^ for n > — 1 if initial values 
satisfy 

a;o,x_2 7^x_i. (34) 

Relative to the multiplicative group of all nonzero real numbers, the second order 
equation fl33l) is HDl with form symmetry 

H{u,v) = — , h{s) = -. 
V s 

11 



Making the substitution (!29|) gives the type-(l,l) order reduction 

_ *o 

Sn+l — O-Sn-i Sq — 

Now using (!27j) and ( l28l) we obtain the following formula for solutions of (132 
subject to (IMl) : 

I J- \^ i 7(7+i)/2 ''0 Xo — a; 



a^n = a^o + to > Sott"" ' " , sq 



Example 3. Consider the following variant of Eg. (1321) : 

x„+i = x„ H , a 7^ (35) 

Xn—\ Xn—2 

subject to (1341) . As in Example 2, the HDl form x„ — x„_i gives the SC factorization 

tn+l = ; to = Xq — X_i, t_i = X_i — ,T_2 (36) 

tn-1 

Xn+1 = ^n+1 + Xn- (37) 



Unlike ( l33l) . Eg. (1361) is not HDl. But a straightforward calculation shows that 
every solution of (1361) has period 6 as follows 

ato a? 0? at_i 1 , , 

t-i,to,-—,- — ,—,^—(- (38) 

t-l t„i To to J 

Thus we may use f l371) and fl28l) to calculate the corresponding solution of fl35l) 
explicitly: For each n, there are integers 5„ > and < p„ < 5 such that n = 65„+p„. 
If cr is the sum of the six numbers in (l38l) then the explicit solution of Eg. (1351) may 
be stated as 

Xn = Xo + Y.l=lti = XO + -^{n - pn) + Er=n-p„+l^i 

where every U is in the set (l38|) and the last term is zero if p„ = 0. 



12 



Remark. ( The full triangular factorization property) 

The equation in Example 2 has an interesting extra feature: it can be fully SC 
factored as a system of first order difference equations 

This system is triangular in the sense that each equation is independent of the 
variables in the equations below it. For a general discussion of the periodic solutions 
of systems of this type see [2] , [TT] . 

If a difference equation of order k + 1 has the property that it can be factored 
completely into a triangular system of first order difference equations then we say 
that the difference equation has the full triangular factorization property or that it 
is FTF. Clearly, every HDl equation of order 2 is FTF but it is by no means clear if 
all HDl equations of order 3 or greater are FTF. For instance, it is not obvious that 
Eq. ( I35l) in Example 3 does in fact have the FTF property. The difficulty there is due 
to the non-HDl nature of ( l36l) which leads to a form symmetry that involves complex 
functions (see Example 5 below). 

For non-HDl equations, the FTF property is not clear even for equations of order 
2. But in Corollary 1 in the next section we show that every linear non-homogeneous 
equation of order A; + 1 is FTF with a complete factorization into a triangular system 
of linear non-homogeneous first order equations. 

5 Separability and type-(l, k) factorization 

The class of HDl functions does not include certain familiar functions. For example, 
the linear non-homogeneous function (j){u,v) = au + bv + c is HDl relative to the 
group of all real numbers under addition only when a + b = 1; it is HDl relative to 
the group (0, oo) under ordinary multiplication only when c = and a,b > 0. These 
restrictions suggest that a proper study of order reducible form symmetries for linear 
difference equations does not belong in the context of HDl equations. 

In this section we define a class of equations that properly includes all linear non- 
homogeneous difference equations with constant coefficients as well as some other 
interesting non-HDl equations. Before discussing this class, recall that a type-(l,A;) 
equation has a SC factorization with factor of order m = 1 and cofactor of order k. 
Therefore, the form symmetry is a scalar function that may be written as 

H{uo, ...,Uk)=uo* h{ui, ...,Uk). 
13 



Now we define a function : 6*^'+^ — i> G to be separable (or algebraically factorable) 
relative to G if there are k + 1 functions (pj : G ^ G, j = 0,1, . . . ,k such that for all 
uo,. . .,Uk e G, 

0(mo, . . . ,Uk) = 4>o{uo) * ■ ■ ■ * 4>k{Uk)- 

Note that every linear non- homogeneous function is trivially separable relative to 
every additive subgroup of the complex numbers C. The rational function 0(m, v) = 
au^/v is separable relative to the group of nonzero real numbers under multiplication 
for every integer p but it is HDl relative to the same group if and only ii p = 2. The 
exponential function (f){u, v) = y^'^-^'^-™ is separable relative to the group of non-zero 
real numbers under multiplication but it is not HDl relative to that group. 

5.1 Additive forms 

We define Eq.(|Hl) to be separable if every function /„ is separable relative to the 
underlying group G. In this section we consider the following separable version of (IHl) 
over the group of complex numbers C under addition: 

Xn+i = a„ + 0o(a;„) + 0i(x„_i) H (t^ki^n^k)- (39) 

with 

x_j,an&C, 0j:C^C, i = Q,l,...,k. (40) 

It is not strictly necessary for the sake of applications that the maps (j)j be defined 
on all of C (indeed, in most applications they are defined on the set M of real numbers) 
but we make a strong assumption to reduce the amount of technical details in this 
article. The use of complex numbers is necessary because form symmetries of ( l39l) 
may be complex even if all quantities in PUI) are real (this happens in particular for 
linear equations). 

The next result from [19] shows that Eq. ( |39l) has an order- reducing form symme- 
try if one of the k + 1 functions (po, . . . ,(f)k can be expressed as a particular linear 
combination of the remaining k functions. The form symmetry in this case gives a 
type-(l, k) order reduction of fl39|) . 

Theorem 3. Assume that there is a constant c G C such that the functions (po, ■ ■ ■ ,4>k 
in Eq.^3^ satisfy 

c''+h - c'^Mz) - c''-^Mz) c(f)k-i{z) - Mz) = for all 2. (41) 

Then / f5^) has the following form symmetry 
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H{zo, zi,..., Zk) = Zo + hi{zi) H h hk{zk) (42) 

where 

h,{z)=c^z-c^-'Mz) 0j-i(^), J = l,---k (43) 

The form symmetry in (^^ and (^^ yields the type-{l, k) order reduction 

Zn+l=an + CZn, Zq = Xq + hi{x ^i) ^ h/lfe(x_fc) (44) 

Xn+1 = Zn+1 - hi{Xn) /lfc(x„_fe+i). (45) 

Remark. Note that the factor equation (H4l) has order 1 and the cofactor ( 145|) has 
order k in this case. For reference, we note that (HTl) and (H3!) imply the following 

c/ifc(^) = (t)k{z). (46) 



A significant feature of Eg. (1451) is that it has the same form as (l39l) . Thus if the 
functions hi, . . . ,hk satisfy the analog of ( HTl) for some constant c' G C then Theorem 
3 can be applied to (H5l) . The next result exploits this feature by applying Theorem 
3 to a linear non- homogeneous equation repeatedly until we are left with a triangular 
system of first order linear equations. 

Corollary 1. The linear non-homogeneous difference equation of order k + 1 with 
constant coefficients 

Xn+l + boXn + biXn-1 H h bkXn-k = "n (47) 

where bo, . . . ,bk, «« ^ C has the FTP property and is equivalent to the following 
triangular system of k + 1 first order linear non-homogeneous equations 

Zo,n+l = Otn + Co2;o,n; 
^1,71+1 = 20,n+l + CiZi^n 



Zk,n+1 ^fc— l,n+l ~r C^Zk^n 

in which Zk^n = Xn is the solution of Eq.^^ and the constants cq, ci, . . . ,Ck are the 
eigenvalues of the homogeneous part of ^1\) , i.e., roots of the characteristic polyno- 
mial 

P(z) = z^+^ + boz^ + biz^-^ + ■■■ + bk^iz + bk. (48) 
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Proof. Defining (pjiz) = —hjZ for j = l,...k and applying Tlieorem 3 above 
yields the SC factorization 

Xn+1 = ^0,n+l ~ l^lflXn — ■ ■ ■ — jSi^k-l^n-k+l (49) 



where Cq satisfies fHT|) 

c^'^^z + c^h^z + c^Q~^hiz H h cohk-iz + hkZ = Q 

for all z G C, i.e. Cq is a root of the characteristic polynomial P in fH8|) . Further, the 
numbers /5i j are given via the function hj in (HSl) and fH6l) as 

Alternatively, the numbers /3ij may be calculated from the recursion 
A,i = Co/?i,j-i + &j, j = 1, . . . fc - 1, /5i,o = Co + 6o, Co/?i,fc-i = -&fc-i- (50) 

Next, since Eq.fHOl). i.e., Xn+i + A.o^^n + " — I" Pi,k-iXn-k+i = ^o,n+i is of the same 
type as ( l47l) . Theorem 3 can be applied to it to yield the SC factorization 

Zl,n+1 = Zo^n+1 + CiZi^n 
^n+l = ^l,n+l ~ P2,0'^n — ■ ■ ■ — P2,k-2Xn-k+2 

in which Ci satisfies (HI!) for (l49l) . i.e., the power is reduced by 1 and coefficients 
adjusted appropriately as in 

Ci + A,oC^' + /3i,ict"2 + . . . + p,f^_^co + (3i,k-i = 0. 

Now we show that Ci is also a root of P in (jUj). Define 

Pi{z) =z'' + /3i,o2;^"' + ■ ■ ■ + /3i,fc-2^ + (3i,k-i 

so that Ci is a root of Pi. If it is shown that {z — cq)Pi{z) = P{z) then Pi divides P 
so ci is a root of P. Direct calculation using (1501) shows 



(z - co)Pi(2;) = z^+^ + /?i,o;z'' + A,iz'^ + ■ ■ ■ + (5i,k-iz 

- C^Z^ - Co/?l,o^^~^ Co/5i,fc_22; - Co/5i,fc-l 

= Z^^^ + (co + hQ)z^ + (CoA,o + feO^'^-' ■ ■ ■ + (CoA,fc_2 + 6fc_2)2: 

- c^z^ - co/3i,o2;^~"^ co/5i,fc-22: - co/5i,fc-i 

= pf^). 
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Therefore, the above process inductively generates the system in the statement of 
this corollary. 

Remarks. (Operator factorization, complementary and particular solutions) 

1. The triangular SC factorization of Corollary 1 is essentially what is obtained 
through operator factorization. If Exn = Xn+i represents the forward shift operator 
then as is well-known, the eigenvalues factor the operator P{E) with P defined by 
(HHjl : i.e., (H7|) can be written as 

{E - Cq){E -Ci)---{E - Ck)Xn-k = Oln- (51) 

Now if we define 

{E -ci)---{E - Ck)xn-k = yo,n (52) 

then ( ISTI) can be written as 

yo.n+l - Coyo,n = Oln 

which is the first equation in the triangular system of Corollary 1 with ?/o,n = ^o.n- 
We may continue in this fashion by applying the same idea to fl52l) : we set 

{E -C2)---{E - Ck)Xn-k = yi,n 

and write fl52|) as ?/i,n+i — Ciyi,n = ^o,n which is the second equation in the triangular 
system if yi „ = Zi^n-i- The reduction in the time index n here is due to the removal 
of one occurrence of E. Proceeding in this fashion, setting i/j^n = Zj^n-j at each step, 
we eventually arrive at 

(E — Ck)Xn~k = yk-l,n =^ Xn+l-k = Z/fc-l,n + CkXn-k- 

Thus, with yk-i,n = -^fc-i.n-fe+i the preceding equation is the same as the last 
equation in the system of Corollary 1. 

2. With Corollary 1 we may obtain the eigenvalues and both the particular solu- 
tion and the solution of the homogeneous part of fHTl) simultaneously without needing 
to guess linearly independent solutions, namely, the complex exponentials. We indi- 
cate how this is done in the second order case k = 1 which is also representative of 
the higher order cases. First, for a given sequence s = {s„} of complex numbers and 
for each c G C, let us define the quantity 



cr„ 
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and note that for sequences s,t and numbers a,b E C, an{as + bt;c) = acr„(s; c) + 
ban(t; c), i.e., an{-, c) is a linear operator on the space of complex sequences for each 
n > 1 and each c G C. Further, if s„ = a6" then it is easy to see that 

Now, if /c = 1 then the semi conjugate factorization of ( 1471) into first order equations 
is 



2;„+i = a„ + co2:„, zo = a;o + (co - bo)x^i (54) 

a;„+i = ^n+i + ciXn- (55) 



A straightforward inductive argument gives the solution of ( 1541) as 

Zn = zoCq + cr„(a; cq). (56) 



Next, insert (1561) into ( 1551) . set 7„ = 2;„+i and repeat the above argument to obtain 
the general solution of (H7|) for k = 1, i.e., a;„ = Xqc" + (J„(7; Ci). If Ci 7^ Cq then from 
( 153|) we obtain after combining some terms and noting that Jq = Zi = a^ + CqZq, 

Xn = ] Cq+ [xo q + an{(J (a; Cq); Ci). 

V Co - Ci / V Co - Ci / 

where a' = {cr„+i(a;; co)}. We recognize the first two terms of the above sum as giving 
the solution of the homogeneous part of (H71) and the last term as giving the particular 
solution. In the case of repeat eigenvalues, i.e., ci = cq again from (l53l) we get 

Xn = [xqCq + (cto + coZo)n]co~^ + cr„(a'(a; co); cq). 

5.2 Multiplicative forms 

As another application of Theorem 3 we consider the following difference equation on 
the positive real line 

yn+l = Pn1po{yn)i'liyn-l) ' " ■^fc(yn-fc), (57) 

(3n, y-j e (0, oo), ^pj : (0, oo) ^ (0, oo), j = 0, . . . k. 



Taking the logarthiin of Eq. (l57|) changes its multiphcative form to an additive 
one. Specifically by defining 



J/r; 



Inyn, Vn = e""", ln/?„ = a„, 0j(r) = ln?/'j(e''), j = 0,...k, re 



we can transform (1571) into (139!) . Then Theorem 3 implies the following generalization 
of the main result of 1 181. 



Corollary 2. Eq.(52) has a form symmetry 

H{to, ti, . . . , tk) = tohiih) ■ ■ ■ hk{tk) (58) 

if there is c e C such that the following is true for all t > 0, 

MtfMtf"---Mt) = t'''^'- (59) 

The functions hj in ^5^ are given as 

/i,(t)=r>o(t)-^"'---^,-i(t), J = l,...k (60) 

and the form symmetry in l[58\) and (E3) yields the type-{l, k) order reduction 

rn+i = /5„<, ro = yohiy-i) ■ ■ ■ hk^y-k) (61a) 

hi[yn> ■ ■ ■ hk[yn-k+i> 

Example 4 (A simple equation with complicated multistable solutions). Equations of 
type fl39l) or fl57|) are capable of exhibiting complex behavior, including the generation 
of coexisting stable solutions of many different types that range from periodic to 
chaotic. As a specific example of such multistable equations consider the following 
second-order equation 

Xn+i = Xn-ie"-''--^"-', a;_i, xq > 0. (62) 

Note that Eq. flB21) has up to two isolated fixed points. One is the origin which is 
repelling if a > (eigenvalues of linearization are ±e"/^) and the other fixed point 
is X = a/2. If a > 4 then x is unstable and non-hyperbolic because the eigenvalues 
of the linearization of (!62|) are —1 and 1 — a/2. The computer-generated diagram in 
Figure 1 shows the variety of stable periodic and non-periodic solutions that occur 
with a = 4.6 and one initial value x_i = 2.3 fixed and the other initial value Xq 
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Figure 1: Bifurcations of solutions of Eq. ()62p with a changing initial value; a = 4.6 is fixed. 



changing from 2.3 to 4.8; i.e., approaching (or moving away from) the fixed point x 
on a straight hne segment in the plane. 

In Figure 1, for every grid value of xq in the range 2.3-4.8, the last 200 (of 300) 
points of the solution {a;„} are plotted vertically. In this figure, stable solutions with 
periods 2, 4, 8, 12 and 16 can be easily identified. All of the solutions that appear in 
Figure 1 represent coexisting stable orbits of Eq. fl^^ . There are also periodic and non- 
periodic solutions which do not appear in Figure 1 because they are unstable (e.g., 
the fixed point x = 2.3). Additional bifurcations of both periodic and non-periodic 
solutions occur outside the range 2.3-4.8 which are not shown in Figure 1. 

Understanding the behavior for solutions of Eq. 062p is made easier when we look 
at its SC factorization given by (16 lap and (161bp . Here k = 1 and 



Mt) 



-t 



iJi{t) = te \ Pn = e" for all n. 



Thus fl59l) takes the form 



e-'He-' 



f foralH>0 
e' for alH > 
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The last equality is true if c = —1, which leads to the form symmetry 



-1 , U\-l 1 . rr/ N Wo 



h{t) = r'Mty =—-,^ H{uo,ui, 

te~* Uie 



and SC factorization 



rn+i = — , ro = Xohi{x^i) = ^^-— (63) 

^"+1 = , T \ = rn+iXne~''". (64) 



All solutions of (1^ with vq ^ e"/^ are periodic with period 2: 



e" 1 f ,To a;_ie" ^- 



^0, — 



ro I I a;_ie -^-^ xq 



Hence the orbit of each nontrivial solution {x„} of (!62|) in the plane is restricted 
to the pair of curves 

e" 
^i(t) = — te"* and ^2(t) = rote"*. (65) 

ro 

Now, if x_i is fixed and xq changes, then tq changes proportionately to xq. These 
changes in initial values are reflected as changes in 'parameters in flMl) . The orbits of 
the one dimensional map 6te~* where 6 = tq or e"/ro exhibit a variety of behaviors 
as the parameter h changes according to well-known rules such as the fundamental 
ordering of cycles and the occurrence of chaotic behavior with the appearnce of period- 
3 orbits when h is large enough; see, e.g., [3], [5], [14j, ^22j. Eq. flM]) splits these 
behaviors evenly over the pair of curves (jHSD as the initial value Xq changes; see Figure 
2 which shows the orbits of (!62l) for two different initial values xq with a = 4.6; the 
first 100 points of each orbit are discarded in these images so as to highlight the 
asymptotic behavior of each orbit. The splitting over the pair of curves ^i,.^2 also 
explains why odd periods do not appear in Figure 1. 

Example 5. We now combine different types of form symmetry to show that Eq. (!35!) 
in Example 3 has the FTF property. If we assume that 

a > 0, Xo > x_i > x_2, Xo > (66) 

then the multiplicative group (0, oo) of positive real numbers is the invariant set of 
Eq.fl36l). Since fl36|l is obviously separable over (0, cxd), we check equality fl59|) in 
Corollary 2 with ipo{t) = t and ipi{t) = 1/t: 

f-^ = f' for all t > 0. 
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Figure 2: Two of the orbits in Figure 1 shown here on their loci of two curves .^i, .^2 in the 
state space. 

This condition holds if c^ — c + 1 = 0. The quadratic has complex roots 



c± 



l±iV3 



so by Corollary 2, Eq. fl5^ has a form symmetry 

H{uo,Ui) = Uohi{ui), where hi{t) = t^+t^^ = t~ 
and an SC factorization 

rn+i = ar^+, ro = tohi{t_i) 

tn+l = ^n+ltn ■ 



(67) 
(68) 



The three equations (jSZ]), (IB51) and ^S7\ establish that Eq. flS^ has the FTF prop- 
erty with a factorization 



r„ 



ar, 



c+ 



C++C- 



tn+1 — CLT^ t^ , 

It is noteworthy that in spite of the occurrence of complex exponents, this system 
generates positive solutions from positive initial values. This fact may seem less 
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surprising if Eq. (136!) is transformed into a linear equation (with complex eigenvalues) 
by taking logarithms as in the beginning of this section. 

Remark. Under the added restrictions (l66l) . Corollary 2 may also be applied to 
Eq. ( 1331) of Example 2 to obtain an SC factorization using the separable type of form 
symmetry. Is this SC factorization different from that in Example 2? To see that 
they are in fact the same, let ipo{t) = t^ and 'ipiit) = 1/t in the equality fl59|) and 
require that 

This holds if c is a root of the quadratic c^ — 2c + 1 = 0, i.e., c = 1. Thus using 
(!60|) we calculate the form symmetry as 

H{uo,ui) = uohi{ui), where hi{t) = t^t"'^ = t^^. 

This is just the HDl form symmetry giving the same SC factorization as in Ex- 
ample 2. 
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